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Abstract 

Avermectins are potent and popular veterinary pharmaceuticals used globally to 
fight parasites of livestock and humans. By disturbing ion channel transport 
through the membrane, avermectins are effective against endo- and ectoparasitic 
round and horsehair worms (Nematoida), insects, or ticks (Arthropoda), but not 
against Plathelminthes, including flatworms (Trematoda) and tapeworms 
(Cestoda), or segmented worms (Annelida). Unfortunately, excreted avermectins 
have strong nontarget effects on beneficial arthropods such as the insect commu- 
nity decomposing livestock dung, ultimately impeding this important ecosystem 
function to the extent that regulators mandate standardized eco-toxicological 
tests of dung organisms worldwide. We show that the ancient phylogenetic 
pattern and qualitative mechanism of avermectin sensitivity is conserved and 
compatible with most recent phylogenomic hypotheses grouping the Nematoida 
with the Arthropoda as Ecdysozoa (molting animals). At the species level, we 
demonstrate phylogenetic clustering in ivermectin sensitivities of 23 species of 
sepsid dung flies (Diptera: Sepsidae). This clustered 500-fold quantitative varia- 
tion in sensitivity may indicate recent lineage-specific responses to selection, but 
more likely reflects pre-existing genetic variation with pleiotropic effects on eco- 
toxicological responses to pollutants. Regardless, our results question the com- 
mon practice in eco-toxicology of choosing single test species to infer detrimental 
effects on entire species communities, which should ideally assess a representative 
taxonomic sample. 



Introduction 

Avermectins are a potent and popular class of veterinary 
pharmaceuticals used globally to fight parasites of live- 
stock and even humans (Omura 2008). By disturbing 
ion channel transport through the membrane, avermec- 
tins are effective against endo- and ectoparasitic round 
and horsehair worms (Nematoida), insects, or ticks 
(Arthropoda), but not against Plathelminthes, including 
flatworms (Trematoda) and tapeworms (Cestoda), or 
segmented worms (Annelida) (Omura 2008). However, 
excreted avermectins also have strong nontarget effects 



on beneficial arthropods such as the insect community 
decomposing livestock dung, ultimately impeding this 
important ecosystem function (Wall and Strong 1987; 
Floate et al. 2005; Jochmann et al. 2011; Lumaret et al. 
2012). Systematic disturbance of dung decomposing 
organisms, primarily beetles, flies and earthworms, by 
toxic anthropogenic veterinary pharmaceuticals excreted 
by livestock has produced sufficiently negative effects 
that regulators mandate environmental risk assessments 
worldwide (VICH 2004; EC 2009). Eco-toxicological lab- 
oratory tests follow strict standardized guidelines (OECD 
2008, 2010) and involve few, often single test species 
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such as, for dung dwellers, the yellow dung fly Scathoph- 
aga stercoraria (Diptera: Scathophagidae), the face fly 
Musca autumnalis (Diptera: Muscidae), or the dung bee- 
tle Aphodius constans (Coleoptera: Scarabaeidae; Rombke 
et al. 2009, 2010a; OECD 2008, 2010). The implicit 
assumption in eco-toxicology is that a single test species, 
typically of temperate origin, represents and typifies the 
sensitivities of various organisms to a particular toxic 
substance present in natural communities (EFSA 2010; 
Rombke et al. 2010b). This ignores that even closely 
related species may have very different natural sensitivi- 
ties to toxins. 

Recent studies of natural populations (e.g., McKenzie 
and Clarke 1988; French-Constant et al. 1993; Weill 
et al. 2004) or experimental evolution (e.g., Lopes et al. 
2008; Vogwill et al. 2012) have shown that resistance to 
frequently man-made toxins can evolve and spread very 
fast by natural selection, mediated either by point muta- 
tions within single genes or by exploiting pre-existing 
genetic variation in physiological pathways with some- 
times unexpected pleiotropic effects on toxin sensitivity. 
In the simplest scenario, organisms may be naturally 
resistant to particular toxins (e.g., antibiotics) simply 
because they already bear a particular gene or genetic 
mechanism that may or may not have another function, 
and which happens to become co-opted or exapted 
(Gould and Vrba 1982) by chance to affect sensitivity to 
a novel drug (Weill et al. 2004; Allen et al. 2010; Walsh 
and Duffy 2013). For example, studies of Drosophila not 
only revealed that DDT resistance involves over-expres- 
sion of the Rst(2)DDT locus, but also that fruit flies 
bearing Rst(2)DDT alleles display pre-existing cross- 
resistance to two other drugs used for pest control in 
agriculture (e.g., Daborn et al. 2001). 

Comparative phylogenetic analyses can provide valuable 
insights about the relative importance of recent episodes of 
selection versus pre-existing natural genetic variation 
underlying drug resistance. Recent or ongoing episodes of 
strong selection are expected to diminish the signal of 
shared common ancestry, which in the extreme will annihi- 
late any correlation between phylogeny and the trait (Losos 
1999). Moreover, when comparing the sensitivity of organ- 
isms to particular toxins, we would expect a pattern in 
which populations or species inhabiting the same geo- 
graphic region share a higher propensity for being drug 
resistant, even if they are not closely related to each other, 
due to their common history of exposure to a particular 
toxin. In contrast, strong phylogenetic clustering of toxin 
sensitivity would rather suggest co-option of pre-existing 
genetic variation prior to exposure to the toxin. 

We have previously assessed the sensitivity, in terms 
of mortality, to the prominent and widely used parasit- 
icide ivermectin for 23 temperate (Europe, North 



America) and tropical (Asia, Central America) sepsid 
dung fly (sub)species in an ecotoxicological study 
(Blanckenhorn et al. 2013). All species used are ecologi- 
cally bound to fresh feces of large mammals or, rarely, 
waterfowl, in which larval development takes place. In 
nature, these species can be collected from fresh cow 
dung, often several species on the same cowpat, which 
is also the breeding substrate of choice in our labora- 
tory. Here, we map our ecotoxicological data onto the 
most recent sepsid phylogeny based on 10 genes (Lei 
et al. 2013) to explore the compatibility of the ancient 
phylogenetic pattern of avermectin sensitivity with most 
recent phylogenomic hypotheses (Burda et al. 2008; 
Telford and Copley 2011) grouping the Nematoida 
together with the Arthropoda as Ecdysozoa (molting 
animals). We further investigate more recent phyloge- 
netic patterns of ivermectin susceptibility among closely 
related sepsid flies. Our results suggest that ivermectin 
sensitivity in sepsid flies evolved largely independently 
of exposure to avermectins. 

Methods 

We worked with 23 sepsid species that were identified 
using Pont and Meier (2002) and SepsidNet (Ang et al. 
2013). The sample included differentiated North Ameri- 
can and European subspecies (populations) of Sepsis neo- 
cynipsea and S. punctum (Puniamoorthy et al. 2012). For 
7 temperate Sepsis species (Sepsis cynipsea, S. fulgens, 
S. neocynipsea, S. orthocnemis, S. punctum, S. thoracica, 
and S. violacea; Table SI), we collected multiple (2-5) 
populations, which originally served for estimating 
within-species variation in sensitivity (Blanckenhorn et al. 
2013); these data were averaged in our main analysis 
here. Depending on availability, we also collected two or 
more of our test species from the same sites (i.e., popu- 
lations): Vienna, Austria (5 spp); Zurich, Switzerland (6 
spp); Sierra Nevada, Spain (2 spp); Tartu, Estonia (2 
spp); Umbria, Italy (3 spp); Calabria, Italy (2 spp). This 
allows comparing systematic geographic variation in iver- 
mectin sensitivity across species, which may reflect local 
adaptation to past drug exposure. All flies used were 
originally caught by ourselves in the wild at the various 
sites specified in Table SI, on or around cow dung. Ide- 
ally, multiple (>10) females for each species were caught 
at each site of origin, which were subsequently kept in 
our laboratories in Zurich and/or Singapore for multiple 
generations using standard methods (e.g., Puniamoorthy 
et al. 2012), either in large population containers or in 
smaller containers as iso-female lines (i.e., offspring of 
one field-caught female). 

Laboratory tests were performed in three temporal 
blocks (2008, 2009, and 2011) using standard toxicological 
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methods with six ivermectin concentrations (0.21, 0.66, 
2.08, 6.57, 20.75, and 65.7 /ig ivermectin/kg dung fresh 
weight) plus water and acetone controls (five replicates per 
concentration; complete methods and toxicological results 
are reported in Blanckenhorn et al. 2013; Table SI). Sepa- 
rately for each species or population, lethal ivermectin con- 
centrations causing 50% larva-to-adult mortality (LC50) 
plus their 95% confidence intervals were estimated by 
regressing logit-transformed emergence proportions 
against log 10 (ivermectin concentration) (ToxRat 2003). For 
comparison, we additionally estimated adult body size 
(head width in mm) and egg-to-adult development time 
(in days) of all flies in the control (i.e., water and acetone) 
treatments (see Table SI), two life history traits that are 
generally subject to natural selection (e.g., Blanckenhorn 
et al. 2007), and which are nonlethally reduced and pro- 
longed, respectively, by ivermectin (Rombke et al. 2009; 
Blanckenhorn et al. 2013). 

We used Mesquite (Maddison and Maddison 2011) to 
map mean ivermectin LC50, body size, and development 
time data as continuous traits onto the most recent sepsid 
phylogeny (Lei et al. 2013) to reveal potential phylogenetic 
clustering of ivermectin sensitivity at the (sub)species taxo- 
nomic level. Phylogenetic clustering of our three traits was 
independently tested via 1000-fold randomization of trait 
values without replacement across species on the given 
phylogeny, yielding a P-value indicating how extreme the 
actual distribution of species on the phylogeny (Fig. 2) is 
relative to all other possible, randomly assembled distribu- 
tions. To investigate potential concerted evolution by natu- 
ral selection of life history traits with ivermectin sensitivity, 
we further computed correlations between the three traits 
using independent contrasts (CAIC: Purvis and Rambaut 
1995). 

At the ancient phylogenetic level, we mapped avermectin 
sensitivity onto the most recent and now widely accepted 
phylogenomic hypothesis grouping the Nematoida together 
with the Arthropoda as Ecdysozoa (molting animals: Burda 
et al. 2008; Telford and Copley 2011). 

Results 

The modern sister group relationship of Arthropoda and 
Nematoida (Ecdysozoa hypothesis) suggests a single origin 
of avermectin sensitivity in animals, whereas the formerly 
hypothesized sister group relationship of Arthropoda and 
Annelida (Articulata hypothesis) would require assumption 
of multiple origins (Fig. 1; adapted from Burda et al. 
2008). 

At the species level, we found up to 500-fold variation 
in the mean LC50 among the tested related sepsid flies 
(Fig. 2). The variation in sensitivity is nonrandomly dis- 
tributed on the sepsid phylogeny as revealed by trait 
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Figure 1 Traditional (top) and new (bottom) phylogenetic hypotheses 
for the ancient relationships between invertebrate taxa (adapted from 
Burda et al. 2008). The pattern of ivermectin sensitivity (red = sensitive 
and green = nonsensitive bars) more parsimoniously agrees with the 
most recent phylogenomic hypothesis. 

value randomization without replacement (P < 0.002). 
At least two clades, Archisepsis plus Microsepsis from 
Central America and the species of the Sepsis punctum 
group that are widespread in Eurasia and North Amer- 
ica, have significantly lower ivermectin sensitivities, as 
indicated by higher LC50 values in Fig. 2. Body size is 
also nonrandomly distributed (i.e., clustered: P < 0.032), 
while development time is seemingly randomly distrib- 
uted on the phylogeny (P < 0.154; Fig. 2). Importantly, 
neither life history trait correlates with ivermectin 
sensitivity (|r| < 0.3; P > 0.2), so their patterns evolved 
independently. 

When comparing ivermectin sensitivity (logit-transformed 
LC50 value) among sites while controlling for species in a 
two-way anova for the subset of seven species col- 
lected from more than one site in Europe (seven sites), as 
well as when considering all species, no systematic geo- 
graphic variation was apparent F 6jl4 = 1.38, P = 0.287 (site 
by species interaction removed due to nonsignificance; the 
same results when including all species and sites: 
F 12i i4 = 1.57, P = 0.208). For the subset, variance compo- 
nent analysis revealed approximately 4.5 times more vari- 
ance among species (4.8) than was observed among sites 
(1.1). 
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Saltella sphondylii (0.20 19.89 1.10) 
Archisepsis armata {6.20 17.15 1.20) 
Archisepsis diversiformis (1.92 20.15 1.14) 
Microsepsis mitts (24.31 14.14 0.77) 
Microsepsis armillata (24.58 15.32 0.83) 
Meroplius fukuharai (0.14 17.29 0.98) 
Dicranosepsis emiliae (0.24 17.55 0.87) 
Sepsis dissimilis (0.11 17.65 0.75) 
Sepsis lateralis (0.80 1 6.53 1 .07) 
Sepsis thoracica (0.36 12.50 1.09) 
Sepsis cynipsea (0.31 1 3.60 1 .03) 
Sepsis neocynipsea EU (0.23 17.00 1.13) 
Sepsis neocynipsea NA(1.15 14.40 1.12) 
Sepsis flavimana (0.05 21.82 0.95) 
Sepsis duplicata (0.09 20.82 0.79) 
Sepsis secunda (1 .33 22.58 0.83) 
Sepsis violacea (0.89 21.00 0.97) 
Sepsis fulgens (2.59 16.70 0.99) 
Sepsis orthocnemis (5.49 18.20 0.97) 
Sepsis latiforceps (11.44 14.81 1.07) 
Sepsis punctum EU (5.76 19.10 1.24) 
Sepsis punctum MA (11. 40 13.90 1.06) 
Themira minor (1 .261 1.33 0.90) 



Figure 2 Ivermectin sensitivity (lethal concentration at which 50% of 
individuals die LCB0, in fig ivermectin/kg fresh dung), mean develop- 
ment time (days), and mean head width (mm) (in parentheses) for 23 
sepsid dung fly (sub)species mapped onto their phylogeny (females and 
males combined; cf. Table S1). Tree branches are colored according to 
the LCB0 values (see legend) to visualize the nonrandom sensitivity pat- 
tern. 



Discussion 

Our study shows that eco-toxicological responses to pol- 
lutants, exemplified by the LC50 values used here, can 
and should be treated as traits subject to evolution that 
may follow phylogenetic patterns. At the ancient phyloge- 
netic level, the pattern of avermectin sensitivity of inver- 
tebrates conforms to the modern phylogenomic 
hypothesis grouping the Nematoida together with the 
Arthropoda as Ecdysozoa (molting animals), rather than 
the traditional hypothesis grouping Annelida and Arthro- 
poda as Articulata. Retrospectively, phylogenetic relation- 
ships can thus explain why avermectins unfortunately do 
not work against major parasitic helminth groups such 
as flatworms and tapeworms (Plathelminthes), but only 
against roundworms (Nematoida), and why they at the 
same time negatively affect Arthropoda, which comprise 
numerous nontarget beneficial species such as the beetles 
and flies decomposing livestock dung treated here: Only 
the two latter groups molt, and apparently avermectins 
disturb ion channel transport through the membrane 
especially during molting (Omura 2008). Fortunately, 
avermectins also do not strongly affect annelids such as 
earthworms (Lumaret et al. 2012), which also play a very 



important role as decomposers of vertebrate dung 
(Holter 1979; Lumaret et al. 2012), and this is again con- 
sistent with the most recent Ecdysozoa hypothesis 
depicted in Fig. 1. These patterns could not have been 
predicted by the previously accepted Articulata hypothe- 
sis, showing that phylogenetic information is useful for 
predicting large-scale patterns of resistance to toxins. 
Eco-toxicology thus supports recent findings in phyloge- 
nomics (Telford and Copley 2011), and important new 
insights were gained here from analyzing resistance data 
in a phylogenetic context. 

A highly significant phylogenetic signal in ivermectin 
sensitivity was also revealed in a clade of 23 dung-dwell- 
ing sepsid flies by trait randomization among species on 
our phylogeny. Two subclades independently show lower 
susceptibility to ivermectin (Fig. 2), the Neotropical 
Microsepsis and Archisepsis species, as well as several 
Eurasian Sepsis species of the punctum group (S. violacea, 
fulgens, orthocnemis, latiforceps, punctum: Pont and Meier 
2002). As outlined in the Introduction, this phylogenetic 
pattern can be explained either by natural selection 
against ivermectin sensitivity, that is, by greater adapt- 
ability, or by pre-existing genetic variation that evolved 
in a different context sometime in the distant past and 
which later pleiotropically affected toxin resistance in 
some groups, but not others (i.e., an exaptation: Gould 
and Vrba 1982). Ivermectin has been in use only for 
about 40 years (Omura 2008), corresponding to roughly 
500 sepsid fly generations with a generation time of ca. 
3 weeks. That insects with large effective populations can 
rapidly evolve resistance to insecticides is a well-known 
phenomenon from pest control (McKenzie and Clarke 
1988; French-Constant et al. 1993; Weill et al. 2004; 
Lopes et al. 2008; Vogwill et al. 2012), and this could 
also be the case in sepsid flies. Indeed, a substantial frac- 
tion of the total variation in ivermectin resistance in sep- 
sid flies can be attributed to population differentiation 
(Blanckenhorn et al. 2013). However, there are several 
reasons why this scenario seems unlikely. The less sensi- 
tive species in the punctum group regularly co-occur on 
the same pastures with several other Sepsis species that 
are very sensitive to ivermectin (e.g., S. cynipsea or flavi- 
mana; Fig. 2). Recurrent substance treatment in such 
regions should simultaneously select against sensitivity of 
all sympatric species from various lineages, not only in 
some. Explaining the obtained pattern by natural selec- 
tion for ivermectin resistance would require the parallel 
fixation of mutations in several reproductively isolated 
species or populations within very short time, producing 
geographic patterns of sensitivity. Yet we found that var- 
iation between species in ivermectin sensitivity is about 
4-5 times greater than variation between sites commonly 
harboring populations of several of these species, 
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although we readily admit that our taxon and popula- 
tion sampling was opportunistic and not ideal for sys- 
tematically examining geographic variation in 
susceptibility. Unfortunately, we have no information 
about ivermectin use at the various sites, as the sub- 
stance is not controlled everywhere. We know only that 
ivermectin use varies locally, even by farmer, and that 
the flies' dispersal distance likely transcends and homog- 
enizes any regional patterns of substance use. Therefore, 
while the nonrandom phylogenetic sensitivity pattern we 
found does not principally exclude recent rapid evolu- 
tion of resistance to man-made toxins by natural selec- 
tion, we do not think that it is the most parsimonious 
explanation for our data. 

There are additional reasons why a phylogenetic expla- 
nation of the obtained pattern is more reasonable. It is 
known that organisms can be 'naturally' resistant to 
novel toxins without ever having been exposed to them 
(Weill et al. 2004; Allen et al. 2010; Walsh and Duffy 
2013). Thus, it is probable that the nonrandom pattern 
of ivermectin sensitivity in sepsid dung flies evolved as a 
by-product correlated with some unknown trait that 
changed in the distant past, exemplifying a preadapta- 
tion or exaptation (Gould and Vrba 1982) conferring a 
selective advantage only post hoc in pastoral habitats cre- 
ated by humans. The connection to the ancient phyloge- 
netic level of ivermectin sensitivity to molting (Fig. 1), 
an evolutionary novelty at that time, can perhaps also be 
interpreted as such an exaptation, in this case conferring 
a disadvantage post hoc. Some further, circumstantial but 
anecdotal evidence supports our evolutionary interpreta- 
tion. First, the species most likely exposed to ivermectin 
should be the cow dung specialists Sepsis cynipsea and 
S. neocynipsea that are most abundant in Europe and 
North America, respectively, which however are two of 
the most sensitive species both in the laboratory 
(Blanckenhorn et al. 2013; Fig. 2) and in the field world- 
wide (Madsen et al. 1990; Floate 1998; Iwasa et al. 
2005). Second, one of the lesser susceptible Archisepsis 
species (Fig. 1; Table SI) was collected in a nature 
reserve on an island in Panama (Barro Colorado) that 
has never been used agriculturally. Third, ivermectin has 
been isolated from the producing fungus Streptomyces 
avermectinus in Japan about 40 years ago (Omura 2008), 
so the common ancestors of the two least sensitive sep- 
sid groups in Fig. 2 could have plausibly evolved iver- 
mectin resistance in the distant past in response to 
natural contact with that fungus. This is again unlikely, 
however, because the substance, and the fungus, so far 
has been found only this one time and nowhere else 
(Omura 2008). 

In summary, we conclude that the nonrandom pattern 
of ivermectin sensitivity only affecting molting animals 
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(Ecdysozoa in Fig. 1), as well as that found for sepsid 
dung flies (Fig. 2), likely indicates exaptations (Gould 
and Vrba 1982) at two different levels: first conferring a 
selective disadvantage to all sensitive molting animals, 
and in a second step, a selective advantage to some spe- 
cies post hoc in pastoral habitats created by humans. Sus- 
ceptibility to natural or man-made toxins therefore need 
not necessarily be adaptations, as is often assumed (e.g., 
Tack et al. 2012). Furthermore, at a more practical level 
in the realm of toxicology, the strong phylogenetic signal 
of ivermectin sensitivity we found implies that any par- 
ticular species cannot possibly be representative when 
chosen for assessing toxicity of substances affecting the 
dung community. This makes choice of any test species 
in the context of environmental risk assessments of vet- 
erinary pharmaceuticals particularly delicate (Blancken- 
horn et al. 2013). One obvious solution to the problem 
is to mandate use of several test species based on diverse 
phylogenetic taxon sampling, or even the dung commu- 
nity as a whole (Floate et al. 2005; Jochmann et al. 
2011), thus extending the registration process already 
required for these drugs (VICH 2004). In this context, 
the issue of regionalization of substance control should 
also be considered, as is already being discussed (EFSA 
2010). 
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Additional Supporting Information may be found in the online version 
of this article: 

Table SI. Proportion of flies emerged from the control treatments 
(water and acetone) for all 47 tests using 21 sepsid species, estimated 
lethal concentration at which 50% of the flies died (LC50, in /(g ivermec- 
tin/kg wet dung) with 95% confidence limits, and mean ± SE develop- 
ment time and body size for male and female sepsids (Nm, Nf = total 
number of males, females) {species averages with SD in italics). 
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